Enhanced Tensor-Force Contribution in Collision Dynamics 
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The tensor and spin-orbit forces contribute essentially to the formation of the spin mean field, 
and give rise to the same dynamical effect, namely spin polarization. In this paper, based on time- 
dependent density functional calculations, we show that the tensor force, which usually acts like 
a small correction to the spin-orbit force, becomes more important in heavy-ion reactions and the 
effect increases with the mass of the system. 

PACS numbers: 25.70. Jj, 21.60.Jz, 21.30.Fe 
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INTRODUCTION 



The tensor force, which is necessary to explain the 
properties of the deuteron, attracts special attention re- 
cently, because it has turned out to play an essential 
role in the existence limit of exotic nuclei, as well as 
the nuclear shell structure far from the /3-stability line 
(for example, see [3-S|). An important feature is that 
the spatial average of the tensor force is exactly equal to 
zero, so that its effect is spatially localized. On the other 
hand, the spin-orbit force, which is necessary to explain 
the large spin polarizations of scattered nucleons, plays 
a crucial role in the nuclear shell structure. The origin of 
the tensor force can be found in the one-pion exchange 
potential, and that of the spin-orbit force in the relativis- 
tic aspects of quantum dynamics. 

Thus the tensor and spin-orbit forces are quite differ- 
ent in their origins, while resulting in the same dynam- 
ical effect, namely, spin polarization. Spin polarization, 
which arises mostly from the spin-orbit force, sponta- 
neously takes place in the early stage of heavy-ion reac- 
tions, and affects the equilibration process to a large ex- 
tent. As long as the microscopic time-dependent Skyrme 
energy density functional (Skyrme-EDF) calculations are 
concerned, the appearance of spontaneous spin polariza- 
tion even in central collisions between /3-stable nuclei was 
shown, and its origin was clarified to be the time-odd part 
of the spin-orbit force [ij]. Therefore, the enhancement or 
reduction of spin polarization gives an ideal framework 
to pin down the properties of the tensor force in collision 
situations. 

In this paper, the role of the tensor force in heavy-ion 
reactions is investigated based on time-dependent den- 
sity functional calculations with explicitly implemented 
tensor force, where the time-odd part of the spin-orbit 
force is also fully taken into account. Special attention is 
paid to the effect of the tensor force on time evolution. 
As a result, some information on the importance of the 
contribution from the tensor force in heavy-ion reactions 
is presented. 



THEORETICAL FRAMEWORK 
Mean field due to spin-orbit and tensor forces 

The contribution of the tensor force, whose role was 
underestimated and mostly neglected for a long time, was 
substantially studied in the context of Skyrme-EDF only 
recently £s|— lS| - Here we begin with the functional form 
of the tensor and spin-orbit forces in Skyrme-EDF. Let 
p, er and J represent the number density, spin density, 
and spin-orbit density, respectively. The contribution of 
the tensor and spin-orbit forces to the energy density 
functional has the form 



W,(r) • (-i)(V x „) 



(1) 



where q = n,p (n and p stand for neutron and proton, 
respectively). W g (r), which is called the form factor of 
the spin mean field, is decomposed into the contributions 
from spin-orbit and tensor forces. 



W 9 (r)=W 9 "(r)+W^ (r), 



(2) 



where Wj s (r) and W^(r) denote the form factors of 
spin-orbit and tensor mean fields, respectively. The con- 
tribution of the spin-orbit force to the functional is 
represented by 

W£ s (r) = |W (Vp(r) + Vp,(r)) + f(ti - t 2 )3 q {v). 

Note that the second term on the right-hand side, whose 
contribution in collision situations was discussed in [TO ], 
has never been taken into account for some modern 
Skyrme parameterizations such as SLy4d and SKM*, be- 
cause it makes fitting spin-orbit splittings more difficult. 
Although there are several versions of the tensor force, 
we are concerned with the natural tensor force only. Its 
contribution to the energy functional is represented by 



WT(r) = aJ g (r) + /3.V(r) 



(3) 



with q' = n,p satisfying q ^ q' , according to Stancu- 
Brink-Flocard [l2| , where the unique contribution of the 
tensor force can be found in J g '(r) to W q (r). The full 
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FIG. 1: (color online) Time evolution of 40 Ca + 40 Ca at the bombarding energy 130 MeV (c.m.). Snapshots of the density are 



shown in a fixed square (40 x 
force used is SkM* + SV-tls. 



40 fm ) on the reaction plane, where contour lines are plotted for multiples of 0.04 fm . The 



introduction of the tensor force requires to refit addi- 
tionally the corresponding central-force parameters. Al- 
though the full introduction brings about largely different 
and complicated contributions depending on the choice 
of force parameter sets [l3| , it has been shown to mostly 
result in weakening the natural contribution |14|. Several 
versions of the tensor force are compared in (14j . Here 
we restrict discussion to the tensor force as defined by 
Eq. because the aim is not a discussion of the exis- 
tence limit of exotic nuclei, but rather the general fea- 
tures of the tensor contribution in reactions. It is read- 
ily seen that the effect of the tensor force corresponds 
to a quantitative modification of that due to the spin- 
orbit force. Accordingly, the contribution of the tensor 
force should be discussed in association with the spin- 
orbit force. 

Tensor-force contribution in collision situations 

A framework for measuring the effects of the tensor 
force is presented with a focus on collision dynamics. 
Concerning the spin polarization, it is reasonable to be- 
gin with a discussion of spin-orbit coupling. It is defined 
by the scalar triple product 

L • S = -ih (r x p) ■ (er + cr') 

(4) 

= — ih r • (p x (er + er')) , 

where er and er' denote the spins of the two nucleons. In 
collision situations r x p is related to the impact parame- 
ter. Comparing Eqs. ([T]) and (j4|), W g (r) in Eq. (JTJ) plays 
the role of the vector r in Eq. Q , where the momentum 
p is replaced approximately by V in the Skyrme-EDF. 

In order to evaluate the contribution of the tensor force 
to spontaneous spin polarization, we introduce a proper 
theoretical setting of heavy-ion collisions. Our starting 
point is that the tensor and the spin-orbit forces are lo- 
calized effects, which are not easy to compare in collision 
dynamics, if there is some similarity in their localized 
patterns. Let the reaction plane be (x, z) with the initial 
collision direction z, and the direction perpendicular to 



the reaction plane be y. For simplicity, the spin direction 
of the initial state is assumed to be parallel to the y-axis. 
In this setting, because only the z-component of p and 
the y-component of er are non-zero, we have 

L • S = -iH X [py (a + ct') - p z (cr + cr'), ) 

V 1 (5) 

= ih xp z (cr + u') y . 

We see that only the x-component of the vector r, and 
thus the ^-component of W g (r) play a role. In this set- 
ting, the role of the tensor force in the spin polarization 
can be evaluated by the corresponding x-component of 
W^(r). Accordingly, the tensor and spin-orbit forces 
can be compared, if there is a certain similarity between 
the x-components of W^(r) and (r) (otherwise at- 
traction or repulsion happen irregularly from place to 
place). Note that their similarity, which will be shown to 
be true, is not trivial. In the following the x-components 
of W^(r) and (r) are simply represented by W^(r) 
and W^ s (r), if there is no ambiguity. 

ROLE OF THE TENSOR FORCE 
Spontaneous spin polarization 

A systematic three-dimensional time-dependent den- 
sity functional calculation is carried out in a spatial box 
48 x 48 x 48 fm 3 with a spatial grid spacing of 0.8 fm, in 
which the Skyrme-force parameter set SV-tls [3| is used 
for the tensor part, and SkM* and SLy4d [l6| for the re- 
mainder including the spin-orbit force: a = 71.102 [MeV 
fm" 5 ] and j3 = 35.142 [MeV fm" 5 ]. The parameter set 
SV-tls was lately introduced in the context of the refit- 
ted tensor force; it is one of the most reliable parameter 
sets in terms of reproducing the contribution of the form 
factor W^(r) of the tensor mean field. The relative ve- 
locity in the collisions is set to 10 % of the speed of light, 
and the initial distance of the colliding nuclei to 20.0 fm; 
their initial positions are (0,0,10) and (0,0,-10). In order 
to pay special attention to the mass-dependent general 
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FIG. 2: (color online) Spin distribution (the spin is projected 
onto the y-axis) of a compound nucleus. A snapshot of a 
composite nucleus, which corresponds to the case at time = 
6.0 x 10~ 22 s in Fig. [TJ is shown on the reaction plane. For 
reference, contours of the density distribution are also shown 
(contour = 0.01, 0.06, 0.11 and 0.16 fm -3 ). 



features, we consider central collisions between identical 
N=Z nuclei: 16 + 16 0, 40 Ca + 40 Ca and 56 Ni + 56 Ni. 
The contributions from J q and J q > in Eq. (J3)) are not so 
different for collisions between N=Z nuclei, therefore the 
parameter dependence mostly arises from the sum of a 
and j3. Some features of the tensor force acting on N = Z 
bound nuclei were studied in . 

Figure Q] shows the time evolution of 40 Ca + 40 Ca re- 
sulting in fusion, where the terms associated with the 
tensor force (SV-tls) are explicitly included. The same 
calculation without the spin-orbit force does not achieve 
fusion. Omitting the tensor force while including the 
spin-orbit force shows no notable difference to the den- 
sity evolution with all force terms included. This suggests 
that large dissipation arises from the spin-orbit force, 
while the tensor-force contribution is definitely small. 
The composite nucleus evolves with a continuing oscil- 
lation; the two nuclei get into contact around time = 4.2 
x 10~ 22 s, and the first full-overlap is achieved at 5.6 x 
10- 22 s. 

Let us consider the y-projection of spin for each single 
nucleon. The spin distribution of the colliding nuclei is 
calculated by their superposition: 

P(t,T) = P (t,Th - P (t,T) U 

where p(t, r)f and p(t, r)j, denote the densities of spin-up 
and spin-down components, respectively. In this defi- 
nition, the density plays the role of weight. The value 
of P(t, x) is positive if the spin-up component is more 
abundant, zero for saturated spins, and negative other- 
wise. As is seen from the presence of (er + a') in Eq. 
the problem of comparing the different role of tensor and 
spin-orbit forces becomes meaningless if spontaneous spin 
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FIG. 3: (color online) Snapshots of the x-component of W 9 (r) 
projected on the reaction plane, corresponding to the case 
at time = 6.0 x 10~ 22 s in Fig. [T] The values are plotted 
separately for the tensor and spin-orbit forces, and for protons 
(q = p) and neutrons (q = n), respectively. 



polarization is absent. Spin polarization appears for all 
the reactions and all the force parameter sets used; e.g., 
in Fig. [5J the presence of spin polarization is shown for 
40 Ca + 40 Ca. As a result, the concept of examining the 
role of the tensor force in the presence of spin polarization 
is valid and will be carried out in the following. 

Figure [2] shows that strong spin polarization is located 
on the edge of the density distribution. The localized 
pattern of the spin structure is complicated, leading to 
a complicated localization of attraction and repulsion 
due to the tensor force. The spin distribution is point- 
symmetric with respect to the origin, which reflects the 
symmetry of the central collision. Note that the spatial 
average of spin polarization for the spin-saturated system 
is equal to zero. 

Comparison between tensor and spin-orbit forces 

Let us begin with the effect of the tensor force in a 
compound nucleus formed briefly after the full-overlap 
situation (time = 6.0 x 10~ 22 s). In case of 40 Ca + 
40 Ca, Fig. [3] compares the x-components of W g (r) for 
the tensor and spin-orbit forces. Both distributions are 
antisymmetric with respect to the z-axis, and have sim- 
ilar distributions but different signs and amplitudes. It 
is clearly seen that the tensor-force contribution is op- 
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FIG. 4: (color online) Snapshots of the x-components of the 
form factors of spin mean field Wj(r) (upper ones) and 
Wp S (r) (lower ones) at time = 6.0 x 10 -22 s are shown in a 
square (30 x 20 fm 2 ) on the reaction plane. The maximum 
amplitude A of the function is shown in the lower right-hand 
side of each plot. 



posite to the spin-orbit force contribution, and amounts 
to less than 10 percent of latter. It follows that the to- 
tal contribution from tensor and spin-orbit force is not 
so different from the contribution of the spin-orbit force 
alone. No significant difference is noticed between the 
values for protons and neutrons. This is expected for a 
collision between N = Z nuclei. 

This difference in sign and the smallness of the ten- 
sor force contribution compared to the spin-orbit force 
contribution is found to hold regardless of the choice of 
force parameter set and the mass of the colliding nuclei. 
This difference in sign, however, did not appear for the 
y and z-componcnts. On the other hand, comparing the 
x-components of WT(r) and W^ s (r) at time = 6.0 x 
10~ 22 s for 16 + 1% 0, 40 Ca + 40 Ca and 56 Ni + 56 Ni 
as shown in Fig. 2J there is a highly noticeable increase 
with mass for the tensor-force contributions, while it is 
only modest for those of the spin-orbit force. 

Let us move on to the time-dependent features of the 
tensor force. For the reaction shown in Fig [TJ the time 
evolution of the ratio between tensor and spin-orbit con- 
tributions 



W^Wf(<) = 



max r (W^,r)) 
max r (W£ s (t,r)) 



(6) 



is shown in Fig. [5] In addition, the corresponding x- 
componcnts of Wj(r) and (r) are also shown at 

times 1.5 x 10" 22 s, 6 x 10" 22 s and 15 x 10~ 22 s. The 
isoscalar dipole mode shown in Fig. [5] suggest that the 
full-overlap is achieved at time = 5.5 x 10 -22 s, and the 
maximal elongation of the composite nucleus at time = 
7.25 x 10~ 22 s. The relaxation of the tensor contribution 
is not strongly correlated with that of the isoscalar dipole 
oscillation (density oscillation towards the fused system) . 
The contribution of the tensor force is quite small before 
the contact time (4.2 x 10~ 22 s), increases after the con- 
tact time, achieves local-maximum at times 6.75 x 10~ 22 



TABLE I 

40 



Enhancement of the tensor-force contribution for 



40 Ca + 40 Ca. Values of Eq. © are calculated for different 
force parameter sets and isospins. 



Parameter set 


Protons (q = p) 


Neutrons (q — n) 


SkM* + SV-tls 


5.94 


6.23 


SLy4d + SV-tls 


6.54 


7.11 



s and 9.00 x 10 - 22 s, and relaxes afterwards. 

Several points should be remarked here. First, the 
tensor-force contribution is enhanced in collision situa- 
tions, being up to 10 times larger than before the con- 
tact time. Second, the opposite sign and the smallness of 
the tensor compared to the spin-orbit contributions are 
apparent during the heavy-ion collision but not before 
contact. The opposite sign means that the contribution 
of the tensor force continues to weaken the spin polariza- 
tion during the reaction. Third, the similarity between 
protons and neutrons is confirmed throughout the reac- 
tion. 

The tensor-force contribution is compared for different 
force parameter sets in Table HI where the enhancement 
is calculated by the ratio 



W^/W£ 5 (i = 6.5 x 10" 22 s) 
WJ/W£ s (i = 1.5 x 10" 22 s) 



(7) 



where W^/W^ s (i) is calculated as shown in Eq. (|6|). 
This table shows that the enhancement is true indepen- 
dent of the choice of force parameter sets, and no signif- 
icant difference exists between protons and neutrons. 



Mass dependence 

As already mentioned, the opposite sign and the small- 
ness of the tensor compared to the spin-orbit force con- 
tribution is valid independent of the mass. Note that a 
calculation using SLy4d + SV-tls showed the same fea- 
tures, hinting that this is probably not strongly force- 
dependent. 

Figure [6] shows the mass dependence of the ratio of 
tensor to spin-orbit contributions for protons and neu- 
tron (Eq. ©), where the values at time = 6.0 x 10 -22 s 
are chosen to calculate the ratio. In all cases, this time 
corresponds to the time briefly after the first full overlap 
and shows a relatively large tensor contribution close to 
the first maximum, so that it is legitimate to compare the 
magnitude for the three cases. While the values are not 
exactly the same for the two parameter sets, they show 
the same trend; the tensor force contribution becomes 
larger for reactions involving a heavier nucleus. For the 
heavier cases, 20 percent contribution from the tensor 
force compared to the spin-orbit contribution is noticed 
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FIG. 5: (color online) The time evolution of the ratio of contributions from tensor force to those of the spin-orbit force is shown 
for protons and neutrons, respectively (upper panel), where the calculated points (at multiples 0.75 x 10 -22 s) are connected 
by 3rd-order spline functions. This reaction corresponds to the case shown in Fig[T] For reference, the time evolution of the 
isoscalar dipole (is-dipole) mode is shown by a dotted line (upper panel). The corresponding reaction-plane snapshots of the 
s-components of W^(r) and Wj s (r) are shown in a square (30 x 15 fm 2 ) (lower panel), where the maximum amplitude A of 
the function is shown in the lower right hand side of each plot. 



TABLE II: Mass dependence of the growth of spin polariza- 
tion (for an explanation see text). For both parameter sets, 
values are normalized by the values obtained for 16 O + 16 O. 



Parameter set 




4U Ca + 4l, Ca 


bb Ni + bb Ni 


SkM* + SV-tls 


1.000 


0.396 


0.260 


SLy4d + SV-tls 


1.000 


0.571 


0.085 



154.390 [MeV fm" 5 ] and (3 = 139.910 [MeV fm- 5 ]). The 
major difference is that its amplitude is actually smaller 
than the spin-orbit force contribution, but reaches as 
much as 50% of the spin-orbit contribution in 56 Ni + 
56 Ni. The difference between the two parameters can be 
related to the largeness of the a and (3 values proposed 
in the Refs. [12,1171 compared to Ref. [14 1. 



(SkM* + SV-tls). This is not a negligible effect consid- 
ering the remarkable spin-orbit splitting in the ground 
states of heavy nuclei. This should have a certain impact 
on superheavy synthesis; the tensor force is suggested 
to play a considerable role in whether a heavy compos- 
ite nucleus is formed successfully or not. On the other 
hand, the spin polarization becomes smaller for reactions 
involving heavier nuclei. The statistical ratio of spin po- 
larization 

max r (p t (t,r) - p_i(t,r)) 



Er(Pt(*.*)+P4.M) 

between time = 6.0 x 10~ 22 s and 1.5 x 10 -22 s, which 
corresponds to the amplitude of spin polarization due to 
the collision, is summarized in Table ITT1 Thus the tensor- 
force contribution tends to survive for the heavier cases, 
while the spin-orbit force contribution decreases sharply 
with mass. Note that there is no serious discrepancy 
between neutrons and protons visible in Fig. [6] 

Finally, the validity of the obtained results is also con- 
firmed by additionally examining an old tensor force pa- 
rameter set proposed by Stancu-Sprung (l2l . 17| (a = 



CONCLUSION 

Based on time-dependent density functional calcula- 
tions with explicitly implemented tensor force, the role 
of the tensor force has been studied in the context of col- 
lision dynamics. It is remarkable that the contribution 
from the tensor force is enhanced in collision situations. 
Its contribution is mass-dependent and has considerable 
influence on reactions involving a heavier nucleus. 

As long as heavy-ion reactions between N = Z identi- 
cal nuclei are concerned, the opposite sign and the small- 
ness of the tensor force contribution compared to that 
of the spin-orbit force has been confirmed independent 
of mass. In particular, the opposite sign means that the 
spin polarization, thus the large dissipation due to the 
spin-orbit force, is reduced by the tensor force. We con- 
clude that the tensor-force contribution is rather impor- 
tant in heavy-ion reactions with respect to the magnitude 
of dissipation. The results presented in this paper give 
a solid starting point for future researches clarifying the 
role of the tensor force in heavy-ion reactions involving 
exotic nuclei, where the drastically different contribution 
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FIG. 6: (color online) The ratios between tensor and spin- 
orbit force contributions for protons (left panel; q — p) and 
neutrons (right panel; q — n) as functions of the mass of the 
composite nucleus. The values at time = 6.0 x 10 -22 s are 
chosen to calculate the ratio. 



from J q and J q i in Eq. ([3]) might play a significant role. 
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